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ABSTRACT 

The  challenges  of  waste  management  and  increasing  fuel  crisis  can  be  addressed  simultaneously  with  the 
production  of  fuel  from  plastic  wastes.  Global  warming  and  waste  management  policies  have  forced  for  the  use  of  alternative 
fuels  on  engines.  The  production  of  fuel  from  plastic  wastes  will  indeed  tackle  the  environmental  pollution  problem  of  waste 
plastic  management  in  the  land  fills.  Plastics  being  derived  from  petrochemical  source  has  higher  amount  of  hydrocarbon 
which  yield  oil  with  high  calorific  value.  Engine  tests  have  been  carried  out  using  neat  Waste  Plastic  Oil  and  blends  of  Waste 
Plastic  Oil  in  proportions  of  25%,  50%,  and  75%  with  diesel  as  fuel.  Combinatorial  Mathematics  Based  Approach  (CMBA) 
has  been  adapted  to  choose  the  optimum  blend  for  superior  performance  of  the  engine.  Engine  performance  combustion  and 
emission  characteristics  under  different  loads  for  different  test  fuels  are  discussed.  The  results  of  CMBA  and  experimental 
tests  showed  that  WP025  is  the  optimum  blend. 
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INTRODUCTION 

With  ever  increasing  dependence  on  fossil  fuels,  its  depletion  is  zipping  and  zooming  day  by  day.  The 
shortage  of  fossil  fuels  can  be  overcome  with  utilization  of  alternative  fuels  on  diesel  engines.  Waste  plastic  oil  is  one 
among  the  various  alternative  fuels  that  can  be  used  either  as  neat  fuel  or  as  diesel  blend.  The  raw  material  cost  for 
production  of  oil  from  plastic  wastes  is  zero  as  they  are  available  at  free  of  cost  [1],  Hence,  waste  plastic  oil  (WPO)  is 
acclaimed  to  be  a  promising  alternative  fuel.  Many  researchers  have  proclaimed  that  the  properties  of  waste  plastic 
oil  are  close  to  diesel  [2,  3].  For  all  biodiesels,  the  viscosity  is  higher  than  that  of  diesel  even  though  they  are 
esterified  [4,  5].  But,  WPO  has  lower  viscosity  compared  to  other  biodiesels  and  slightly  higher  than  that  of  diesel. 
Hence,  no  need  of  esterification,  also,  little  bit  of  higher  viscosity  improves  the  engine  lubricity  for  extending  the  life 
of  engine. 

Multi  Attribute  Decision  Making  (MADM)  is  a  discrete  approach  to  identify  the  optimum  alternative, 
where,  predetermined  and  limited  numbers  of  alternatives  are  present.  Various  MADM  methods  are  Simple  Additive 
Weighting  (SAW),  Weighted  Product  Method  (WPM),  Analytical  Hierarchy  Process  (AHP),  TOPSIS,  Taguchi,  Grey 
Relational  Analysis  (GRA),  Graph  Theory  Matrix  Approach  (GTMA)  etc.  Jadhav  [6]  investigated  the  engine 
performance  with  Mangifera  Indica  Methyl  Ester  and  adapted  Taguchi  grey  relational  analysis  to  optimize  BSFC  and 
NOx.  Six  input  parameters  with  five  levels  were  considered  to  form  L25  orthogonal  array.  Fuel  injection  nozzle 
geometry  was  found  to  be  the  most  influencing  parameter.  Razmjooei  and  Nagarajan  [7]  adapted  Response  surface 
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method  to  optimize  the  operating  parameters  to  simultaneously  reduce  NOx  and  particulate  matter.  Dimitriou  et  al  [8] 
studied  the  effects  of  injection  rate  with  EGR  combinations  and  bio-diesel  fuel  rates.  Taguchi  orthogonal  array  with  one 
way  analysis  of  variance  have  been  used  to  estimate  the  importance  of  chosen  factors  and  to  find  the  optimum  parameters 
for  reducing  the  emissions  without  loss  reducing  the  power.  Kadam  and  Jadhav  [9]  adapted  Taguchi  and  multiple 
regression  analysis  to  optimize  vibration  and  performance  on  a  diesel  engine.  L16  array  was  formed  with  Taguchi 
technique  and  multi  regression  analysis  was  used  to  find  the  best  relationship  between  the  input  and  output  parameters. 

Earlier  studies  on  WPO  and  its  blends  as  alternative  fuel  established  the  virtues  of  WPO  in  reducing  the  exhaust 
gas  emissions.  Likewise,  a  lot  of  research  has  been  conducted  in  the  past  to  choose  optimum  alternative  utilizing  GRA, 
RSM,  Taguchi  approach  and  so  on.  A  definite  study  or  a  tool  addressing  the  selection  of  optimum  blend  which  is  simple, 
systematic  and  logical  scientific  method  is  needed.  The  present  study  envisages  analysis  of  engine  performance  with 
various  blends  of  WPO.  Combinatorial  Mathematical  Based  Approach  is  adapted  to  select  the  most  desired  WPO  blend 
that  can  be  suitably  tuned  to  obtain  superior  performance  out  of  the  engine. 

2.  MATERIALS  AND  METHODS 

2.1.  Engine  Setup  and  Procedure 

The  engine  tests  were  conducted  on  single  cylinder  (3.7  kW)  constant  speed  (1500  rpm)  four  stroke  diesel  engine. 
An  eddy  current  dynamometer  was  coupled  to  the  engine.  Various  instruments  like  AVL  H12D  miniature  pressure 
transducer,  AVL  364  angle  encoder,  AVL  Digas  444  gas  analyzer,  smoke  meter  of  AVL  make  were  used  to  measure  the 
desired  parameters.  AVL  INDI  MICRA-602-T10602A  Version  V2.5  software  was  used  for  online  analysis  of  engine 
performance.  The  ignition  delay  period,  starting  of  combustion,  heat  release  rates  and  peak  cylinder  pressure  were  obtained 
from  Pressure  Vs  Crank  angle  diagram.  The  technical  specifications  of  the  engine  are  shown  in  Table  1.  Tests  were 
conducted  at  five  load  conditions  i.e.  3A,  6A,  9A,  12A  and  15A  which  correspond  to  20%,  40%,  60%,  80%  and  100%  of 
load.  Engine  was  run  for  10  min  at  every  load  to  stabilize  and  the  readings  were  recorded.  For  each  refueling  with  other 
blends,  engine  was  run  for  30  min  to  consume  the  fuel  entrapped  in  the  fuel  lines.  Repeatability  was  ensured  by  replicating 
the  experiments  thrice.  The  overall  approach  is  given  in  Figure  1 . 


Figure  1:  Engine  Setup  and  CMBA  Methodology 
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Table  1:  Engine  Specifications 


Make  And  Model 

Kirloskar,  AV 1 

Number  of  Cylinders 

1 

Bore,  mm 

80 

Stroke,  mm 

110 

Compression  Ratio 

16.5:1 

Rated  Power,  Kw 

3.7 

Rated  Speed,  Rpm 

1500 

Fuel  Injection  Timing,  ubtdc 

23 

Fuel  Injection  Pressure,  Bar 

215 

Loading  Device 

Eddy  current  dynamometer 

3.  COMBINATORIAL  MATHEMATICAL  BASED  APPROACH 


Combinatorial  Mathematical  Based  Approach  (CMBA)  is  a  Multi  Attribute  Decision  Making  method,  which  is  an 
integration  of  combinatorial  mathematics  matrix  function  and  Analytical  hierarchy  process.  The  objective  of  CMBA 
methodology  is  to  identify  the  factors  that  influence  the  engine  operating  parameters  and  get  the  most  suitable 
combination. 


The  stepwise  procedure  of  CMBA  is  explained  as  follows, 


Step  1:  Decision  Matrix 


The  decision  matrix  consists  of  alternatives,  attributes  and  relative  importance  [10],  In  this  study.  Load  and  WPO 
blends  are  considered  as  alternatives.  BSFC,  BTE,  NOx,  HC,  CO,  Smoke  and  EGT  are  the  attributes.  Eq.  (1)  shows  the 
decision  matrix  in  general  form. 


M  = 


c  c  c 

'-'11  •••  ••• 


IB 


c„ 


c 


iB 


c  c  c 

Wl  •••  Aj  •••  AB 


(1) 


Where,  Cy  is  the  performance  value  and  A,  B  are  the  number  of  alternatives  and  attributes.  The  attributes  are 
either  subjective  or  objective. 

Step  2:  Normalization 

Normalization  makes  all  values  of  attributes  to  be  non-dimensional.  The  normalization  procedure  adopted  in  this 
study  is  as  follows, 


maX  ;  (Cij  ) 

min,(c,) 


;  when  the  j“  attribute  is  found  to  be  beneficial 


;  when  the  j“  attribute  is  found  to  be  non  beneficial 


Eq.  (4)  shows  the  decision  matrix  with  normalized  values. 


(2) 

(3) 
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In  the  present  study,  BTE  is  beneficial,  whereas  BSFC,  N0X,  HC,  CO,  Smoke  and  EGT  are  non  beneficial.  The 
normalized  values  of  the  experimental  results  are  shown  in  Table  2. 


Table  2:  Normalized  Values  of  Experimental  Results 


Exp 

Fuel 

Load 

BSFC 

BTE 

NOx 

HC 

CO 

Smoke 

EGT 

No. 

(%) 

(kg/h  kW) 

(%) 

(ppm) 

(ppm) 

(%) 

(%) 

(0C) 

1 

WP025 

20 

0.049 

0.024 

1.000 

0.127 

0.922 

1.000 

1.000 

2 

WP025 

40 

0.345 

0.375 

0.884 

0.392 

0.813 

0.914 

0.759 

3 

WP025 

60 

0.633 

0.660 

0.652 

0.608 

0.719 

0.741 

0.578 

4 

WP025 

80 

0.885 

0.869 

0.446 

0.907 

0.422 

0.513 

0.387 

5 

WP025 

100 

1.000 

1.000 

0.251 

1.000 

0.000 

0.216 

0.156 

6 

WPO50 

20 

0.035 

0.009 

0.989 

0.082 

0.953 

0.984 

0.960 

7 

WPO50 

40 

0.336 

0.322 

0.862 

0.302 

0.859 

0.902 

0.709 

8 

WPO50 

60 

0.619 

0.592 

0.636 

0.526 

0.766 

0.690 

0.492 

9 

WPO50 

80 

0.881 

0.796 

0.417 

0.840 

0.484 

0.404 

0.271 

10 

WPO50 

100 

0.982 

0.932 

0.219 

0.981 

0.047 

0.161 

0.070 

11 

WP075 

20 

0.022 

0.004 

0.973 

0.019 

0.969 

0.980 

0.899 

12 

WP075 

40 

0.323 

0.311 

0.826 

0.228 

0.891 

0.894 

0.668 

13 

WP075 

60 

0.611 

0.548 

0.607 

0.485 

0.828 

0.678 

0.447 

14 

WP075 

80 

0.872 

0.751 

0.354 

0.799 

0.516 

0.341 

0.226 

15 

WP075 

100 

0.965 

0.902 

0.148 

0.922 

0.063 

0.090 

0.030 

16 

WPOIOO 

20 

0.000 

0.000 

0.964 

0.000 

1.000 

0.973 

0.854 

17 

WPOIOO 

40 

0.310 

0.285 

0.786 

0.194 

0.922 

0.882 

0.613 

18 

WPOIOO 

60 

0.593 

0.527 

0.571 

0.433 

0.875 

0.667 

0.392 

19 

WPOIOO 

80 

0.850 

0.705 

0.311 

0.765 

0.578 

0.294 

0.166 

20 

WPOIOO 

100 

0.947 

0.869 

0.000 

0.854 

0.188 

0.000 

0.000 

Step  3:  Relative  Importance  Values  for  Attributes 


The  decision  maker  meticulously  analyses  the  attributes  and  assign  relative  importance  values  (off  diagonal 
elements)  using  Analytic  hierarchy  process  scale  [11,  12].  The  pair  wise  comparison  matrix  N,  as  shown  in  Eq.  (5),  is 

1 


formed  with  p^  and  pj;  where  p;j  =  1  when  i=j  and  p  = 


Pn 


1 

Pn 

Pn-  ■ 

■  P  IB 

Pn 

1 

Pn  ■ 

P  2B 

N  = 

Pi  i 

1 

P  IB 

_Pm 

PA2  ■ 

1 

The  scale  for  pair  wise  comparison  is  sh( 

(5) 
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Table  3:  Pair  Wise  Comparison  Scale 


Degree  of  Importance 

Definition 

1 

Equal 

2 

Intermediate  between  1  and  3 

3 

Moderately  preferable 

4 

Intermediate  between  3  and  5 

5 

Strongly  preferable 

6 

Intermediate  between  5  and  7 

7 

Very  strongly  preferable 

8 

Intermediate  between  7  and  9 

9 

Extremely  strongly  preferable 

Pair  wise  comparison  matrix  in  the  present  study  is  shown  in  Eq.  (6), 


BSFC 

BTE 

NOx 

HC 

CO 

Smoke 

EGT 

BSFC 

1 

3 

2 

3 

4 

6 

6 

BTE 

0.333 

1 

2 

2 

3 

3 

5 

NOx 

0.500 

0.500 

1 

2 

3 

3 

4 

HC 

0.333 

0.500 

0.500 

1 

2 

4 

3 

CO 

0.250 

0.333 

0.333 

0.500 

1 

2 

2 

Smoke 

0.167 

0.333 

0.333 

0.250 

0.500 

1 

3 

EGT 

0.200 

0.250 

0.333 

0.500 

0.333 

0.333 

1 

The  consistency  check  [13]  should  be  carried  out  to  ensure  the  judgments  of  relative  importance  values  of 
attributes.  As  per  the  decision  making  method,  consistency  ratio  should  be  less  than  0.1.  For  pair  wise  comparison  matrix 
shown  in  Eq.  (6),  the  consistency  ratio  was  found  to  be  0.061. 

Step  4:  Matrix  for  Attributes  of  Alternatives 

Matrix  for  attributes  of  every  alternative  is  formed  with  normalized  values  as  diagonal  elements  and  off  diagonal 
elements  is  taken  from  Eq.  (6).  This  matrix  is  represented  as  R  in  Eq.  (7). 


dn 

Pxi 

Px3 •  ' 

■■  PXB 

Pix 

d22 

P  23 

P2B 

R  = 

Pm 

d33 

P2B 

Pax 

PA3  • 

d  ab 

In  the  present  study,  the  alternative  selection  attribute  matrix  is  shown  as  Eq.  (8). 
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d\\ 

0.507 

0.308 

0.324 

0.289 

0.310 

0.250 

0.120 

d22 

0.308 

0.216 

0.217 

0.155 

0.210 

0.180 

0.085 

d33 

0.216 

0.217 

0.155 

0.178 

0.120 

0.085 

0.077 

d\\ 

0.145 

0.207 

0.125 

0.090 

0.056 

0.051 

0.054 

^55 

0.103 

0.083 

0.060 

0.056 

0.051 

0.027 

0.036 

d(,(, 

0.125 

0.072 

0.042 

0.051 

0.054 

0.024 

0.017 

d(,e 

Step  5:  Permanent  Function 

The  “index  score”,  is  also  called  as  Permanent  function  (Per  C)  [12]  and  is  calculated  as  follows. 

Per  (C)  = 

n  r  •  x  x . x ^ 

i  =  l  i  =  l  j  =  i+l  A  =  /  + 1 

+  X  X  X  .  X  (PijPjkPu  +  p*PvPji)p,pmp„po-  p,pu 

i  =  l  j  =  i+lk  =  j  +  l  A  =  t  + 1 

A  —  3  A  A  —  1  A  A 

+  [X  X  X  X . X  (Pa p Pn p,k)p,„p„po p,pM  + 

i  =  l  j=i+lk  =  i+ll=i+  2  A  =  /  +  l 

A  -3  A  -1  A  A  A 

X  X  X  X  X  ( PijPjkPuPu  +  Pi,  Plk  Ptj  P  ji)PmPnPo .  P,P„A 

i  =  l  j  =  i  +  \k  =  i  +  \l=  j  +  l  A  =  t  +  \ 

A  -  2  A  -1  A  A  -1  A  A 

HXXXXX  . X  ( P  ij  P  it  P  k,  +  P  it  PtjP  a  X  P  to  P  ,m )  P„  P„ . P,  P,„ 

1  =  1  j=i+lk  =  j  +  l  l  =  1  m  =  l  + 1  A  =  1 +  1 

A-4  A  -1  AAA  A 

+  X  X  X  X  X  .  X  (PijP  jt  Pkl  Pto  Pto  +  Pto  Pto  Pit  Ptj  P  ji)PnPo .  P,  Pm  1 

1  =  1  j  =  i+lk  =  i+ll= i+lm  =  j  +  l  A  =  l  +  1 

A-3  A-l  A  A  A  -1  A  A 

HXXXXXX  .  X  (PijP  jtPttPli  +  P„  Pit  Ptj  P  ji)(  Pmn  Pn„,)Po P  t  P  m 

1  =  1  j  =  i+lk  =  i+ll=  7  +  1  m=  1  n  =  m+ 1  A  =  1  + 1 

A  -5  A-l  A  A  -  2  A-l  A  A 

+  X  X  X  X  X  X  . X  (P  ij  P  jt  P  ti  +  P  it  P  tj  P  a  X  P  to  P  ™  Pm  +  P  In  P  nm  P  nd  )  P„ P,  P,n 

1  =  1  j  =  i+lk  =  7  +  1  1=1  m  =  l  +  l  n  =  m  +  1  A  =  1  + 1 

A  -5  A  A-3  A  A-l  A  A 

+  X  X  X  X  X  X  . X  (  Pij  P  a  )(Pt,  Pit )( P„,„  Pn„,  )pc .  P,P„, 

1  =  1  j=i+lk  =  i+ll=i+2m  =  k  +  l  n  =  k  +  2  A  =  l  +  1 

A-5  A-l  A  A  A  A  A 

+  X  X  X  X  X  X  . X  (PijP  jt  PtiPto  P  mil  Pm  +  Pin  P  nm  Pto  Pit  Ptj  P  ji  )  Po .  P,Pm  ] 

1  =  1  7  =  i  +  l^  =  i  +  l/=i  +  lm  =  i+l  n  =  7'  +  !  A  =  1  +  1 


A  computer  program  was  generated  to  simplify  the  above  calculation. 

Step  6:  Rank  of  Alternatives 

Calculate  the  index  scores  for  all  alternatives  using  step  5  and  sort  either  in  descending  order.  The  alternative  with 
highest  index  score  is  ranked  as  lowest  from  among  the  alternatives. 

4.  RESULT  OF  CMBA 

As  an  example,  index  score  for  alternative  1  is  presented  below. 

The  alternative  selection  attribute  matrix  is  given  as  (Step  4), 
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0.049 

0.507 

0.308 

0.324 

0.289 

0.310 

0.250 

0.120 

0.024 

0.308 

0.216 

0.217 

0.155 

0.210 

0.180 

0.085 

1.000 

0.216 

0.217 

0.155 

0.178 

0.120 

0.085 

0.077 

0.127 

0.145 

0.207 

0.125 

0.090 

0.056 

0.051 

0.054 

0.922 

0.103 

0.083 

0.060 

0.056 

0.051 

0.027 

0.036 

1.000 

0.125 

0.072 

0.042 

0.051 

0.054 

0.024 

0.017 

1.000 

The  index  score  calculated  for  experiment  no.  1  as  shown  in  Eq.  (10)  is  1.2585.  Likewise,  the  index  scores  for  all 
experiments  are  calculated  and  tabulated  to  rank  them  and  are  shown  in  Table  4. 


Table  4:  Rank  of  the  Alternatives 


Rank 

Experiment  No. 

Fuel 

Load (%) 

Index  Score 

1 

5 

WP025 

100 

1.2585 

2 

9 

WPO50 

80 

1.2364 

3 

14 

WP075 

80 

1.2101 

4 

3 

WP025 

60 

1.2029 

5 

19 

WPO  100 

80 

1.1938 

6 

4 

WP025 

80 

1.1821 

7 

8 

WPO50 

60 

1.1749 

8 

10 

WPO50 

100 

1.1660 

9 

13 

WP075 

60 

1.1624 

10 

18 

WPO  100 

60 

1.1464 

11 

15 

WP075 

100 

1.1352 

12 

20 

WPO  100 

100 

1.1127 

13 

2 

WP025 

40 

1.0897 

14 

7 

WPO50 

40 

1.0650 

15 

12 

WP075 

40 

1.0463 

16 

17 

WPO  100 

40 

1.0313 

17 

1 

WP025 

20 

0.9407 

18 

6 

WPO50 

20 

0.9313 

19 

11 

WP075 

20 

0.9178 

20 

16 

WPO  100 

20 

0.9131 

Table  4  shows  that,  index  score  for  experiment  no.  4  is  highest  among  all  experiments.  As  such,  WP025  at  100% 
load  is  found  to  be  the  optimum  blend  for  superior  performance  of  the  engine.  CMBA  methodology  was  adapted  by  [14- 
20]  in  manufacturing  process  and  similar  type  of  findings  were  reported. 

CONCLUSIONS 

The  present  study  is  to  analyze  diesel  engine  performance  combustion  and  emission  characteristics  using  blends 
of  WPO  and  to  find  the  optimum  blend  for  superior  performance  of  the  engine.  The  following  conclusions  are  drawn  out  of 
this  study: 

1.  The  properties  of  all  WPO  blends  are  acceptable  and  favorable  as  fuel  on  diesel  engine. 

2.  Engine  operated  smoothly  not  only  with  WPO  blends  but  also  with  100%  WPO. 

3.  Waste  plastic  oil  blend  of  WP025  at  100%  load  is  found  to  be  the  best  blend  for  superior  performance  of  the 
engine  based  on  performance  and  emission  characteristics.  The  proposed  method  of  CMBA  considers  any  number 
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of  quantitative  and  qualitative  attributes.  Also,  it  offers  objective  and  simple  assessment  approach.  The  concept  of 
permanent  function  provides  complete  information  without  any  loss. 
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